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SUMMARY 

The theory developed by G. B. Whitham (Com- 
munications ov Pure and Applied Mathematics , 
August 1952) for the supersonic flow about bodies 
in uniform flight in a homogeneous medium is 
reviewed and an integral which expresses the (fleet 
of body shape upon the flow parameters in the far 
field is reduced to a form which may be readily 
evaluated for arbitrary body shapes. This ex- 
pression is then used to investigate the e fleet of nose 
angle, fineness ratio , and location of maximum body 
cross section upon the far-fleld pressure jump across 
the bow shock of slender bodies. Curves are pre- 
sented showing the variation of the shock strength 
with each of these parameters. It is found that, for 
a wide variety of shapes having equal fineness ratios, 
the integral has nearly a constant value. Hence, 
to a first order , the pressure jump in the far field is 
independent of the shape and depends only upon the 
fineness ratio. 

INTRODUCTION 

Airplanes operating at supersonic speeds may 
produce shock waves of considerable intensity at 
ground level. The shock pattern travels with the 
airplane, sometimes sweeping over large popu- 
lated areas on the ground below. The atmosphere 
experiences a sudden increase in pressure. as the 
shock wave passes through it. This pressure 
jump across the shock is heard by the observer 
as a sharp explosive-type sound, the so-called 
“sonic boom.” Under some circumstances dam- 
age may result to building components. Thus, 
the sonic boom has become a serious operating 
problem and reliable theoretical prediction of the 
magnitude of sonic booms lias become increasingly 
desirable. A variety of conditions— such . as, 


airplane size, speed, altitude, flight path, winds, 
and atmospheric nonuni formi tics — arc known to 
affect the strength of a sonic boom. In the 
present paper the effects of body geometry are 
investigated . 

One of the more promising theoretical methods 
for predicting the strength of shock waves from 
aircraft at large distances appears to be that 
developed by G. B. Whitham (ref. 1). This 
theory is a rather lengthy modification of the 
linear supersonic theory of slender bodies which 
predicts the location and strength of the shock 
waves emanating from a body in addition to the 
surface pressures and wave drag. Whitham’s 
results, which apply to thickness effects only, may 
be conveniently separated into two parts: the 
prediction of the complicated pressures and shock 
interactions near the body and the greatly simpli- 
fied flow pattern at sufficiently large distances. 
Whitham shows that the calculation of the 
pressure jump across the bow shock in the far 
field requires only a single formula which gives 
the magnitude of the discontinuity in terms of the 
flight Mach number, distance from the flight path, 
and a coefficient which depends upon body shape 
and must be evaluated for each configuration 
under study. 

Whitham’s theory has considerable theoretical 
and experimental verification. When applied to 
two-dimensional problems the method gives a 
first approximation to the result previously 
obtained by K. O. Friedrichs (ref. 2). The 
theoretically predicted far-field strength and 
decay of the bow shock of a supersonic body 
agrees favorably with the experimental evidence 
gathered from studies of bullets in flight (ref. 3), 
wind-tunnel investigations (ref. 4), and flight tests 
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of full-scale airplanes (refs. 5, 6, and 7). On the 
basis of these results it is felt that Whitham’s 
theory will serve as a useful guide in predicting 
the trends of the overpressures associated with 
sonic booms. 

The present paper is concerned with the effect 
of body shape and fineness ratio upon the thick- 
ness-induced pressures in the far field. For back- 
ground the fundamental hypothesis upon which 
Whitham’s work is founded is given, the physical 
basis for the theory is discussed, and then a brief 
development of the equations for the far-field 
conditions is given. The expression for the effect 
of body shape upon the flow parameters in the 
far field is reduced to a convenient form and then 
applied to a number of families of body shapes 
chosen to investigate the effect of fineness ratio, 
location of maximum thickness, and nose angle 
on the pressure jump across the bow shock. The 
reader interested only in the results may proceed 
directly to the application. 

SYMBOLS 

Cb body-shape constant defined in 
equation (24) 


| 


ZttJ 

'o \'y—Z 

™=r 

F(y)dy 

^ max 

largest value which I(y) assumes for 
given body 

l 

length of body 

M 

free-stream Mach number 

n 

integer 

V 

ambient pressure 


free-stream pressure 

> 

II 

1 

00 

R 

radius of body-wake combination 

R max 

maximum radius of body-wake com- 
bination 

R'{ o) 

nose semiangle 

s 

cross-sectional area, i rTl 2 


free-stream velocity 

x,r 

cylindrical coordinates ;%is measured 
along body axis downstream from 
nose, r is radial coordinate meas- 
ured perpendicular to x 

(Xll) max 

distance from nose at which radius 
of body is i? mar 

y&n 

variables 


7 polytropic index, 1.4 

A interval size for numerical integration 

6 local direction of flow 

fj. ' local Mach angle 

<p velocity potential 

Subscripts: 

b body geometry 

w wake geometry 

DESCRIPTIVE OUTLINE OF WHITHAM'S WORK 

WHITHAM’S BASIC HYPOTHESIS 

The theory of reference 1 is a modification of 
the linear supersonic theory of nonlifting slender 
bodies. The basic idea on which the theory is 
based is that the failure of linear theory away from 
the body surface and near shocks is due not so 
much to incorrect prediction of physical quantities 
along the characteristics but to improper placing 
of the characteristic curves in the flow field. More 
precisely, the fundamental assumption upon which 
WhithanVs work is based may be stated as follows. 
Linear theory gives a first approximation to the 
entire flow provided the value it predicts for any 
physical quantity at a given distance from the 
body axis on a straight Mach line through a given 
point on the body surface is interpreted as the 
value at the given distance from the axis on the 
exact characteristic curve pointing downstream 
from the given point on the body surface. Thus, 
the failure of linear theory, in which the charac- 
teristics are straight Mach lines, may be remedied 
by introducing a more exact expression for the 
form of the characteristics which takes their 
bending into account . The curved characteristics 
may converge in various regions of the flow field 
indicating the presence of shocks, and may diverge 
in other regions of the flow indicating expansions. 

THE PRESSURE AS PREDICTED BY LINEAR THEORY 

The physical implications of this basic hypo- 
thesis may be seen with the aid of figures 1 and 2. 
Figure 1 shows the characteristic field and a typical 
pressure trace of a simple slender-body — wake 
combination as predicted by linear theory. The 
wake, shown by the dashed curves, converges 
near the body and rapidly thins to a cylinder of 
approximately constant cross section. The char- 
acteristics — such as, AB, EF, and WZ — issuing 
from various points along the surface and the 
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Figure 1. — Characteristics and pressure as predicted by 
linear theory. 



Figure 2. — Characteristics and pressure as predicted by 
the modified theory. 

characteristics in the upstream uniform flow are a 
series of parallel straight lines, the Mach lines 

x—f3r= Constant, which make an angle ju=sin -1 ^ 

with the body axis. Let I, K, L, . . . R, T 1; T 2 
represent various field points along a line paral- 
lel to the body axis but at some given distance, 
say 7, from it. The atmosphere experiences a 
gradual compression as it passes through the 
forward portion of the characteristic field as- 
sociated with the body (between I and K), and t lie 
ambient pressure begins to rise above free-stream 
pressure. After reaching maximum compression 
at K the atmosphere expands (between K and R) 
and the ambient pressure falls to some minimum 
value at R which is below free-stream pressure. 
Between the maximum and minimum values the 
ambient pressure becomes equal to the free-stream 
pressure at S. After the ambient pressure reaches 
a minimum at R the atmosphere is again com- 


pressed; the ambient pressure slowl}^ rises to free- 
stream pressure. 

THE CHARACTERISTICS AND SHOCKS IN THE MODIFIED 
THEORY 

Conspicuously absent from linear theory are 
shock waves, clearly visible in photographs of pro- 
jectiles in flight, and their associated pressure dis- 
continuities. The first step in introducing the 
si lock waves into the flow consists in deriving a more 
accurate expression for the shape of the character- 
istics and then replacing the Mach lines of linear 
theory by these more correct curves. The shape 
of the characteristics depends upon the local flow 
direction and the local speed of sound. The vari- 
ation of these quantities from their main-stream 
values is neglected in linear theory. In the ap- 
proximation used by Whitham (ref. 1) the char- 
acteristics in the uniform upstream flow remain 
unchanged from the straight Mach lines of linear 
theory, whereas the characteristics issuing from 
the surface of the body and wake are curved, their 
exact shape depending upon the cross-sectional- 
area development of the body. There will, how- 
ever, be one straight characteristic issuing from 
a point near the center of the body. This charac- 
teristic extends to infinity, and everywhere along 
it the ambient pressure is equal to the free-stream 
pressure. The equation of this line is given by 
x—^r—y 0 where the exact value of y 0 depends 
upon the body shape. In general, the curved 
characteristic through a point on the body surface 
is upstream from the Mach line through the same 
point in those regions for which the ambient pres- 
sure is more than the free-stream pressure and is 
downstream from the Mach line in regions where 
the ambient pressure is less than the free-stream 
pressure. 

As a result of this varying shape and curvature, 
the characteristics tend to run together and over- 
lap in some regions of the flow. In such regions 
shock waves arc introduced by the “angle” 
condition: if two regions of supersonic flow are 
separated by a shock then, to the first order in the 
strength, the direction of the shock bisects the 
Mach direction of the two regions of flow. (The 
Mach direction at a point is the outward direction 
making the local Much angle with the local flow 
direction.) The resulting shock waves appear 
in the flow in such a manner as to cut ofl the 
characteristics before they overlap. 
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THE PRESSURE AS PREDICTED BY THE MODIFIED THEORY 

The modifications of the flow field in figure 1 
are shown in figure 2. In this figure the points 
C, H, E, . . . W correspond to the points similarly 
labeled in figure 1. Associated with each of these 
points there is a curved characteristic such as 
CD', EF', UV' which replaces the corresponding 
Mach line of linear theory, CD, EF, and UV in 
figure I. The shocks are shown as the curves 
labeled AD Tv' and XT'X'. The Mach line 
x—pr—y 0 effectively separates the characteristic 
field into two distinct portions, one associated 
with each of the two shock waves. In order to 
obtain the pressure at any point in the field, the 
basic hypothesis stated previously is used. Let 
K', L', M', . . . T' represent field points along a 
line parallel to the body axis and at the given 
distance ? from it. Let the points be chosen so 
that curved characteristics which pass through 
them intersect the body at t lie same places as do 
the Mach lines through the points K, L, M, 

. . . T b T 2 in figure 1 . Then, the basic . lijqjo thesis 
states that the pressure at I/, M', . . . S', X' is 
the same as the pressure predicted by linear theory 
at the points L, M, . . . S, N in figure 1. In 
this manner the pressure may be found at any 
point not on a shock wave. 

Any point on a shock wave is the meeting place 
of two characteristics each carrying a different 
pressure, and hence a sudden increase in pressure 
is experienced as a shock wave is crossed. In 
figure 2 the upstream Mach line through Iv' 
carries with it the free-stream pressure, whereas 
the downstream curved characteristic UK' carries 
the same pressure as point K in figure 1. Thus, a 
pressure jump equal to Ap at Iv is experienced as 
the shock is crossed at K'. Tins method of deter- 
mining the pressure discontinuity along the shocks 
shows that the strength of 1 lie bow shock at first 
increases until maximum strength is reached at 
some distance from the body axis, and thereafter 
decays with distance. Similarly, at T' on the rear 
shock, where the curved characteristics through 
G and W on the body-wake surface meet, the 
excess pressure on the upstream side of the shock 
is the same as that at T L ; the excess pressure down- 
stream is the same as that at T 2 . The pressure 
jump as the shock is crossed at T' is again the 
difference in these two pressures. 

For a body shape which gives rise to a number 
of regions of positive and negative overpressure, 


the characteristic field and the resulting shock 
system will, of course, be far more complicated 
than that described here. However, it is important 
to notice that the modifications of linear theory 
described above tend to smooth out the flow 
pattern between the shocks by "pushing” local 
irregularities away from the Mach line x—(tr=y 0 
and into the shocks as the distance from the body 
increases. Thus, although a system of more than 
two shocks may appear in the neighborhood of 
the body, in general these will coalesce at some 
distance leaving only two shocks extending to 
large distances. In the far field there is an abrupt 
pressure rise across the front shock, a nearly linear 
decline in pressure between the shocks, and a 
rccompression at the rear shock to nearly the free- 
st ream pressure. Tins is the typical N-wavc 
pressure pulse associated with the sonic boom. 

FURTHER ASSUMPTIONS AND RESTRICTIONS TO 
THE THEORY 

In addition to the basic hypothesis discussed, 
the assumptions from which Whitham's theory 
proceeds are the usual requirements of linear 
theory —that is, the body is axisymmotric, slender, 
and pointed at the nose. The upstream field is 
uniform with Mach number sufficiently in excess 
of 1.0 for the bow shock to be attached. The 
flow within the wake is not determined; however, 
the mean boundary of the wake associated with a 
given body is assumed known. The wake is then 
treated as a solid extension of the body; that 
is, there is no flow across the assumed wake 
boundary. The radius of the body-wake combi- 
nation must be continuous although the slope 
may be discontinuous. 

Under these assumptions, reference 1 gives a 
first-oider approximation to the entire pressure 
field due to thickness. A method is described for 
determining the location of the shocks, the pressure 
jump across the shocks, the pressure signature 
between the shocks, and the flow field behind the 
body outside the wake. 

If the projectile is not axisymmotric or has a 
small fineness ratio the shock pattern and pressures 
in the near field may differ appreciably from those 
of a slender symmetric body to which the theory 
is applicable. However, at sufficiently large dis- 
tances from such a body the disturbances will be 
small and the theory may be used to give useful 
information of the far-field conditions. For a 
parabolic body of fineness ratio 5 the trend of the 
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far-ficl(l decay of the bow-shock overpressure as 
predicted by theory agrees well with experimental 
results (ref. 4). 

The assumption that the body is pointed at the 
nose, a usual requirement in linear theory is 
implicit in WhithanTs equations for the perturba- 
tion velocil ies. Erroneous results may be obt ained 
if the theory is applied to a blunt shape. For 
example , the theory is inapplicable to an elliptical 
body (major axis in stream direction) although 
experimental work (ref. 4) shows that the far-field 
pressures are essentially the same as those for a 
parabolic body of the same fineness ratio. 

At low supersonic Mach numbers the bow shock 
will be detached. A region of subsonic flow then 
exists in the vicinity of the body nose. The shock 
in the neighborhood of the nose is then consider- 
ably stronger than in the case of an attached 
shock. It is difficult to justify, on a theoretical 
basis, the use of WhithanTs equation in this 
instance. However, flight test data (ref. 5) ob- 
tained at low supersonic Mach numbers indicate 
that reasonably good agreement is obtained in this 
case. A detached shock also exists ahead of a 
blunt-nosc projectile. There is a real difficulty 
in applying the results to such a shape as already 
noted for the elliptic body. 

THE BASIC EQUATIONS OF THE MODIFIED THEORY 

For the sake of completeness a brief derivation 
is given of the equations for the far-field conditions 
pertinent to the front shock. The numerical 
evaluation of the equation obtained for the pres- 
sure discontinuity will be dealt with in some detail 
subsequently. Certain general results, applicable 
to a variety of axisymmetric problems, arc first 
derived, and then the salient features of the far 
field of a smooth body arc obtained by introducing 
the appropriate approximations for largo values 
of r. 

From the linearized theory of supersonic flow 
it is known (ref. 8) that the form of the perturba- 
tion potential appropriate to a nonlifting body of 
revolution is 


rx-pr 

<p{x 9 r)=— I 


m 


v w-ty-pr* 


At 


where the function /(Q is to be determined from 
the boundary condition of tangential flow at the 
body surface. The perturbation velocities arc 


then obtained by differentiating <p: 


d<p 

.m 

r*-pr 

fit) f/ A 

U dz 

i 

N 

1 

2* 

“V 

Jo V(£“ 

-t) 2 —fi 2 r 2 

^ __d<P x 

m 

1 C X -Z T {x 

d t " 

1 ~dr~~r 

Vz 2 -/3V 2 

r J o 
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For axisymmetric slender smooth bodies, S(x) and 
S'(x) continuous, the condition of tangential flow 
U 

requires that / 0*0 =tt 2 ’ &'(/) where S(x) is the 

Z7T 

cross-sectional area and is the free-stream 
velocity. For small disturbances the body must 
be pointed at the nose, >S r '(0) = 0; hence, /(0)=0 
and only the integral terms need be retained in 
equations (1). 

Following the basic hypothesis outlined in the 
previous section, Whitham replaces x— jSr, the 
characteristic variable of linear theory, by y{x } r). 
The function y(x,r) is determined from the con- 
dition that y(x,r) = Constant he a characteristic 
curve; that is, along the curve 

^=cot(/x+0) (2) 


whore m is the local Mach angle and 6 is the local 
direction of flow. The characteristic variable y is 
defined uniquely by the convention that at any 
point on the body surface y is the value of the 
linear characteristic x—(3r which passes through 
that point. Thus, the equation y=x 0 —(3r 0 defines 
a nonlinear characteristic curve which passes 
through the point x 0} r a on the body surface. For 
bodies of large fineness ratio l/2l\ maXf I3li<^x and 
y^x. Hence, a convenient interpretation of the 
variable y is that it is the horizontal distance 
from the nose at which the extended characteristic 
curve intercepts the body axis. 

The basic results of the modified theory regard- 
ing the pressure at any point in the field, the 
shape of the curved characteristics, and the loca- 
tion of the front shock may be obtained as follows. 
In equations (1) (where now /(0)= 0), if x is 
replaced by y+fir and the integrands arc approx- 
imated for small values of 2///3r, the following 
expressions arc obtained: 


u 


uT 

V2j9r 

► 

V 

,-/3— 

u m * 

p u.J 
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in which 



S"(t) 

Jy-t 


dt 


In the undisturbed flow ahead of the body F(y) = 
0. The pressure at any point in the field is then 
given by equations (3) and the result of linear 
theory as 


, « u yAP F(y) 
7 V m * a/2 


( 5 ) 


The differential equation for the characteristics, 
equation (2), may be approximated to the first 
order in u/U m and vjU m as 




( 0 ) 


The function F{y) defined by equation (4) is 
of fundamental importance in the entire theory. 
Equation (5) shows that, at a given radial distance 
from the body axis, the longitudinal variation in 
pressure is proportional to F{y). The atmosphere 
undergoes compression over those portions of the 
body for which F f (y)^> 0 and expansion where 
F'(y)<C. 0. In general, a shock is formed for each 
value of y for which F"(y) = 0 and F'(y)^> 0 hold 
simultaneously. These shocks usually run to- 
gether at some distance from the body in such a 
manner that only two shocks remain in the far 
field. If the ultimate cross-sectional area of the 
wake Su(<x>) is finite the function F(y) has the 
following properties and physical interpretations: 
(a) The initial compression at the nose: 

F(y) ~2[7?'(0)]V /2 as y-> 0 


where 

y2 j3 m 

equation (6) becomes 


Wlien equation (3) is used, 


clx 

dr 


=fS—kF{y) 


1 

2^ 


( 7 ) 


(b) Recompression outside the wake: 

F{y) ~ y~ zn as y— > <» 

47T 

(c) The balance of compression and expansion: 


If equation (7) is integrated along the character- 
istic 

y=x 0 —fir 0 (S) 

from the point x 0} r o on the body surface to a field 
point £,r, it follows that at large distances on the 
characteristic the following equation must hold: 


F(r})dv = 0 


( 12 ) 


(d) Reexpansion behind a closed body with no 
wake, S w ( oo ) = 0: 

3F 


F{y)~—*rW 


aS y_><S0 


x= fir — IcF iy)r l!2J ry (9) 

Equations (8) and (9) constitute a parametric 
representation of the curved characteristics. The 
two equations may also be interpreted as defining 
the value of y ) and hence of F(y), appropriate to 
any field point x,r , 

With the shape of the characteristics known, 
the “angle” property which must be satisfied 
along a shock is then applied to determine the 
location of the front shock. In terms of the 
parameter y , it is found that along the front 
shock one must have 


x—pr+y—kF(y)r 1,2 

(10) 

'^J/ w 

(11) 


where T r is the volume of the projectile. 

FURTHER SIMPLIFICATIONS APPROPRIATE FOR LARGE 
DISTANCES 

From these general results the simplified equa- 
tions for the conditions at the front shock are ob- 
tained winch are sufficiently accurate in the field 
at large distances from the body-wake combina- 
tion. Tins is the region of greatest concern in the 
prediction of sonic-boom intensities. Since the 
expressions to be obtained are asymptotic results 
for large values of /', they should not be expected to 
give valid results in the neighborhood of the body. 
The conditions under which the results are valid 
arc restated: The body is assumed to be axisym- 
metric, slender, and pointed at the nose (S b (0) = 
jSV( 0) = 0); the projectile-wake surface is smooth 
(S(x) and S'(x) continuous everywhere); and the 
ultimate cross section of t lie wake S w ( °o ) is finite. 
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In addition, it is assumed, as is generally the case, 
that only two shocks reach into the far field. A 
thorough discussion of the structure of the near 
field, the circumstance leading to the appearance 
of three or more shocks in the far field, and the 
methods for dealing with such problems as the in- 
finite slender cone and bodies with discontinuous 
slope, such as the double cone, are given in refer- 
ence 1. These problems require more general 
techniques than are needed for the bod}" shapes 
considered in this paper. 

The far-ficld approximations depend upon the 
fact that between the shocks the only character- 
istics reaching to large distances are those for 
wliicli y~y 0 (fig- 2). The characteristic defined 
by y=y 0 is a straight line and, therefore, from 
equation (9), F(y 0 ) — 0. Figure 2 also indicates 
that in the far field large changes in the radial dis- 
tance along the shocks are equivalent to very small 
changes in y . Hence, with little loss of accuracy, 
equation (11) becomes 

fm-V! w-yl £' FMd ” < 13) 


From equations (13) and (10), the equation for the 
location of the front shock becomes 


x—pr—y<f= 



F(y)dy 


(14) 


From equations (13) and (5), the pressure discon- 
tinuity across the shock is 


A p 

pT 


2 1/4 7 

Vr+I rVi 


vr 


F(y)dy 


(15) 


Equation (15) is the fundamental equation widely 
used in predicting sonic-boom overpressures at 
large distances. 

The pressure at any point between the shocks 
may be determined by approximating equation 
(9) as 


F(y) 


Pr+y 0 —x 

kr 112 


and combining with equation (5) to obtain 

--nw.; «*■-'+») 

Po* (t+1 r 

from which for fixed large values of r the pressure 


between the shock falls linearly at a constant rate 
independent of the body shape. 

NUMERICAL EVALUATION OF THE INTEGRAL 
FOR THE PRESSURE JUMP IN TERMS OF 
CROSS-SECTIONAL-AREA DISTRIBUTION 

It is to be noticed that the detailed behavior of 
the function F(y) is of no concern in the prediction 
of the magnitude and decay of the pressure jump 
across the bow shock in the far field. From equa- 
tion (15) it can be seen that only the value of 

C v ° 

the expression F(y)dy is required. The value 

of this definite integral depends upon shape of the 
body-wake combination under investigation (as 
shown by eq. (4)) and, of course, varies from one 
shape to another. 

The value of the upper limit y 0 is determined 
from the condition that the integral 

I(y)= f F(y)di ? 

Jo 

is a maximum at 7/=7/ a . A necessary condition 
for a maximum at y=y 0 is 

f-JW - 0 

hence, y u is a root of the equation F(y) = 0. How- 
ever, the appropriate root is not necessarily the 
first solution of F(y)= 0 (excluding y 0 = 0) as has 
been indicated in some of the literature concerned 
with WhithanTs method. One must integrate to 
that root which maximizes the integral T{y) in 
order to account for the total strength of any 
subsidiary shocks which may have run together 
near the body to form a single intensified bow 
shock reaching to huge distances. Since the total 
expansion and compression must balance (eq. (12)), 
there is always a root satisfying the required con- 
ditions for body-wake combinations whose ulti- 
mate radius is finite. 

The expression for T(y) will now be reduced to 
a form which avoids determining F(y) explicitly. 
The result to be obtained allows rapid calculation 
of t lie required maximum value and is particularly 
useful in dealing with shapes for which no analytic 
expression is available or for which the analytic 
expression is tedious to handle in an exact manner. 

Write equation (4) in the form 
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Then, the following expression for T(y) is easily 
obtained by integration: 

/('/)=£ F(v)dr,=^ £ S"ft) dt 


Integrating by parts and using the condition 
S"(0) = 0 yields 


m 


p S'(& 

2tt Jo Vy- 




Tliis integral will now be evaluated numerically 
by using an equal-interval Simpson’s rule. Tliis 
rule requires that the upper limit y be an even 
multiple of the interval size which will be denoted 
by A. For the shapes treated in this paper I{y) 
is a well-behaved function which presents no 
difficulty if the integration is carried out analyti- 
cally. However, if t ho integral is to be evaluated 
numerically one must take care in handling the 
singularity of the integrand at the upper limit. 
For this purpose it is convenient to separate the 
integral into two parts 


it may be noted that the last term always involves 
£7(0) which is zero. 

The integral / 2 (y) may similarly be integrated 
by part s : 



The integral in equation (18) is equivalent to the 
expression 

rS(j) + (y-i)S'(y)-S(y) r] . 

J 0 y-S ) 3/2 

which may be simplified by integrating the last 
two terms as follows: 


J- 


at 

3/2 "S 


(y-z ) 3 


+2S'(y)VFH-2-P4 < 19 > 

\V~Z 




i'lu‘1 


r, , S' 

,W= J. n 

w-r si 

J y~ 2 A \y— 


(z) 


dz 


-z 


dz 


If the condition 5(0) =0 is used, hill) may be 
integrated by parts to obtain 

j SQ/—2A) 1 S(Z) 

l0/)_ V2A 2 Jo (y-Z) 3/2 ^ ° 6) 


Applying Simpson’s rule to the integral remaining 
in equation (16) yields 

/,(y)= oM~4^ S(y-2A)~ 3^2 S(y- 3A) 

-4I2 S(y- 4A)- Aj S(y- 5A)- . . .] (17) 

The number of terms in equation (17) depends 
upon the values chosen for y and A; however, 


Thus, subst itut ing equation (19) into equation (18) 
gives 



S(Z)-S(y) 

\¥—z 


-S'(y)\y-z 


I f_5ft) + (y~Z) S' ( y)-S(y) 

-J (y~Z) 3 ' 2 ' 



1 C y S(Z) + (y-Z)S'(y)-S(y) 
zj,-* (y-Z) 3 ' 2 s 


+ S'(vU 2A 


■ 5(y) 
\2A 


S(y-2 A) 
v'2A 


( 20 ) 


Tlie integrand of tlie integral appearing in equa- 
tion (20) is finite over tlie range of integration and 
vanishes at the upper limit. Evaluating this 
integral by using a three-point Simpson’s rule 
gives 

T , , Hy^— 8 rr v ,, > 1 n3^„, OA . 

F{y)=— J2 — VA 5 (?/)— ^-^ [_4~ lS ( y_ 2 A) 

+45(y- A) ~ 13V ^~ 16 5(y)] (21) 

If the results obtained for I\{y) and / 2 (y) (eqs. 
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(17) aiul (21)) arc combined, T(y) is found to be 



- 2§5 S(y- 2A)-^J S(y-3A) 

-^S(y-4A)~ . . .]} (22) 

The number of terms in the series in brackets 
depends upon y and A and terminates with the 
term involving S( 0), which is zero. If no analytic 
expression is available for S(y) } the S'(y) term in 
equation (22) may be eliminated by introducing 
the approximation 


The expression for /(?/) may then be written 


J , s 1 1 r «-nV5 

I{V) - ‘2ir C ^ r A L 4 


- S(y+A )- 


13y'2+lG 


SQj) 


+ 8±HV2 S(3f _ A)+ 2 fi S(y _ 2A) 
S(y- 3A)+ i J 5 S(y- 4A) 


+^»S(y-5A)+ . . . ] (23) 


Equation (23) is the working form desired from 
which the maximum of I(y) maybe readily found. 
It is particularly appropriate for application to 
shapes which cannot be readily expressed analyti- 
cally since it involves only the cross-sectional-area 
distribution. In application a rough calculation 
is first made to isolate the positive peaks in /(?/). 
More than one peak may occur. The peak regions 
arc then investigated more thoroughly to deter- 
mine which is the absolute maximum. The pres- 
sure discontinuity across the bow shock of the 
body is then simply 


Ap_ 2 1/4 t g 1/4 ff 

/W 7 + l r ,/4 A ‘ 


the body length l. If J max, the for- 

miihi for the pressure jump becomes 


Ap p l,i 2 Rj»a; r 


(j)" 


where 




TV’ r 


2 n/4 v'7+f 




(24) 


RESULTS OF APPLICATION TO SPECIFIC BODIES 

In order to investigate the effect of certain body- 
shape parameters upon the magnitude of the bow- 
shock overpressures, the body-shape constant C b 
defined in equation (24) was evaluated for several 
families of body shapes by using equation (23). 
It should be remembered that C b depends upon 
the cross-sectional-area distribution and, hence, 
is a function of such local details as nose angle and 
location of maximum thickness. 

The first family of shapes was chosen to investi- 
gate the effect of varying the location of maximum 
thickness. The body shapes in this group for 
which the maximum thickness lies ahead of the 
center are shown in figure 3. By reversing these 
bodies, four shapes were obtained for which the 
maximum thickness lies behind the center. Since 
the nose angle varies as the location of maximum 
thickness is varied, the effect of the two parame- 
ters has not been completely isolated. 

The effect of nose angle was investigated with 
the bullet-sliape bodies of figure 4. The nose 
angles for these shapes vary from about 5° to 20°. 
The maximum thickness is held fixed at the base. 



It is convenient to normalize S{y) to its maximum Figure 3. — Curves shoving the radius of body shapes 
value Smax ) and normalize the interval size A to having various locations of maximum thickness. 
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Figure 4. — Curves showing the radius of body shapes 
having various nose slopes. 


By normalizing those shapes appropriately the 
group of body shapes shown in figure 5 are ob- 
tained. The nose angle is now constant, about 
7°, and the maximum thickness is again fixed at 
the base; however, the fineness ratio varies from 
8 to 28. 

The shapes shown in figure 6 are similar to the 
shapes obtained at high Mach numbers by apply- 
ing the supersonic area rule to actual airplane 
configurations. The forward portion of the shape 
is the same for all bodies and the location of the 
rearward hump is essentially fixed. Two shocks, 
one from the nose and another from the hump, 
contribute to the strength of the bow shock. The 
five shapes shown in figure 7 arc equivalent bodies 
of revolution for actual airplane configurations. 
These body shapes are included for comparison 
with the idealized shapes treated elsewhere. 

EFFECT OF LOCATION OF MAXIMUM THICKNESS 

Figure 8 shows the variation of the body-sliape 
constant for a family of shapes having various 



Figure 5. — Curves showing the radius of body shapes 
having various fineness ratios. 



Figure 6 . — -Curves shoving the radius of body shapes 
jf similar to those obtained from the supersonic area rule. 


Shope 



Figure 7. — Curves showing the radius of the equivalent 
bodies for several airplane configurations. 



Figure 8. — -Shape coefficient for bodies having various 
locations of maximum thickness. 


locations of maximum thickness. The shapes for 
which the location of the maximum section is 
ahead of xjl~ 1 /2 are determined from the equation 



where 



The constant Ai nor- 
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malizes the right-hand side of the equation so that 

its maximum value on the interval is 

unity. The maximum cross section is located 
{xj V) max =1 2/ (2 -j" 7i) . These shapes are shown in 
figure 3. The shapes for which the maximum 
cross section is behind x/l= 1/2 arc determined 
from the equation 



The location of the maximum section is given by 
(^/0maT=^/(2 + n). The right-hand side of equa- 
tion (26) is also normalized by A\ so that its maxi- 

x 

mum value on the interval O^j^l is unity. The 

shapes given by equation (26) are obtained by 

making the substitution j=l — y in equation (25). 

Thus, for a fixed value of n the shape given by 
equation (26) is identical to that given by equa- 
tion (25), but reversed in the flow. The wake 
assumed in the calculations is defined by R w —0. 
(Although the Von Karman drag is identical for 
two such reversed bodies, it is of interest to note 
that the sonic-boom pressure is different.) 

Figure 8 shows that the body-shape constant 
C b decreases by 20 percent as ( xjJ) max increases from 
0.2 to 0,5. For the reversed shapes C b increases 
as ( x/l) maz increases from 0.5 to 0.8. This increase 
is due to the extreme pointed nose. As ( xjl) ma x 
increases, these shapes behave like bodies having 
a shorter length and a long thin probe pointed 
forward into the oncoming stream. 

Since the nose angle for the nine shapes in this 
family varies with the location of maximum thick- 
ness the effects of these two parameters on the 
value of C b have not been isolated. 


In figure 5, ^ — =4 n. In this case the nose 


angle for all shapes is fixed at 7°, and the fineness 
ratio varies linearly with n. The bod 3 r shapes 
as shown in figures 4 and 5 arc considerably exag- 
gerated in order to display more clearly their area 
development close to the nose. The wake for 
these shapes is taken to be a cylinder extending 
to an infinite distance downstream with a radius 
equal to that of the body base. 

Figures 9 and 10 show a nearly linear increase 
in C b with nose angle and fineness ratio. Although 
these two parameters vary by a factor of 3.5, the 
shape coefficient increases by only 20 percent. 
The increase in C b in figure 9 represents an actual 
increase in the pressure jump due to increasing 
bluntness. However, the pressure jump associ- 
ated with the shapes in figure 5 is proportional to 


2 r„ 


c 

l ^ 6== 4 n an ^ decreases with increasing fine- 
ness ratio. 

For the simple bullet-shape bodies shown in 
figures 9 and 5 the bow shock and the nose shock 



Figure 9. — Shape coefficient for bodies having var'ous 
nose slopes. 


EFFECT OF NOSE ANGLE AND FINENESS RATIO 

The body shapes shown in figures 4 and 5 are 
determined from the equation 


R 


Rm 


= 1 



In figure 4, — =10; this fixes the fineness ratio 

Climax 

of the shapes, but allows the slope at the nose to 
vary according to the relation 



r 'M=To 


Figure 10. — Shape coefficient for bodies having various 
fineness ratios. 
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arc identical. For more complicated configura- 
tions tlic bow shock may be a combination of 
several shocks which coalesce near the body. In 
such cases the nose angle may have even less 
effect upon the bow-shock overpressures than 
indicated in figure 9 as can be seen clearly for the 
next family of shapes discussed where the strength 
of the nose shock is a relatively small contribution 
to the total strength of the bow shock. 


SOME SHAPES SIMILAR TO THOSE OBTAINED BY APPLYING 
THE SUPERSONIC AREA RULE 


The body shapes shown in figure 6 have area 
distributions similar to those obtained at high 
Mach numbers by applying the supersonic area 
rule to actual flight configurations. The shapes 
were obtained from the equation 


(27) 


X 

These shapes are quite flat in the region 0.3^^ 0.5. 


The radius of the body in this region is denoted by 
/iV The important properties of these shapes 
are given in the following table: 


K 

( X}1) max 

R max 


/ I V Ri 1 
\2R *«./ 0r "l *84.4 

Cj, 


Ri 


10 

0. 83 

1. is 

0. 424 

35.8 

1.48 

25 

.87 

1.04 

. 305 

25.7 

2. 17 

40 

.89 

2.15 

. 2325 

19.6 

3. 05 

55 

.89 

2. 00 

.188 

15.9 

3. 95 

70 

.90 

3.18 

.1573 

13.3 

4.83 

85 

.90 

3. 70 

.135 

11.4 

5.70 

100 

.90 

4.22 

. 1185 

10 

0, 57 


As K increases, the location of the maximum cross 
section remains relatively fixed; however, the 
maximum radius increases by a factor of 4 for the 
range of 7v considered here. A given value of K 

Ii 

fixes the value of the ratio jT* . The fineness 


ratio of the shape is then determined from the 


value of the fineness-ratio parameter 



once a value of ^ is chosen. 


The wake used for 


each body in the calculations is shown as the 
broken curve in figure 6. It was chosen so that 

when j-r the flow turns through an angle 

l b4.4 


of 15° after passing over the peak. The flow then 
separates from the rear surface of the projectile 
and asymptotically approaches a cylinder of 
constant cross section. For K=10 the wake is 


simply 7? w — 0. 


When 7T>10, the bow shock which extends into 
the far field from these shapes is a combination of 
two shocks —one from the nose and another from 
the steep incline ahead of the maximum cross 
section — which coalesce at some distance from 
the body. The nose shock has the same strength 
for all shapes; the strength of the second shock 
increases with R rnax . Thus, if y 0 in equation (15) 
were taken as the first root of F{y) = 0 only the 
contribution of the nose shock would be taken into 


account. 


In order to compare the relative strengths of the 
two shocks for a given body, it is convenient to 

introduce the parameter C b —{ ~ /? mqg jC b so that 


A p 
Va> 



For a fixed value of ^ all the 


x 1 

bodies have the same nose shape ahead of and, 

hence, nose shocks of equal strength. The rearward 
hump gives rise to a second shock which augments 
the nose shock in the far field. The differences in 
the size and shape of (lie hump are responsible for 
the variation in C\. The values of C\ are given 
in the preceding table. When the bow 

shock and t lie nose shock are identical. By 
comparing C b for 7^=10 with C b for the other 
shapes it can be seen that the strength of the nose 
shock becomes less important as K increases. 
Thus, any changes in the strength of the nose 
shock due to changes in the nose shape of these 
bodies will be less significant than in those cases 
where the nose shock alone is responsible for the 
fur-field pressure discontinuity. 

Figure 1 1 again indicates only a small variation 
in the coefficient C b in spite of the difference in the 
body shapes. Tn the present case there is a 
20-percent decrease in C b as the fineness ratio 



Figure 11.— Shape coefficient for bodies similar to super- 
sonic-area-rule shapes. 
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varies by a factor of 3.5. This trend in C b is the 
reverse of that shown in figure 10. 

EFFECT OF THE WAKE 

In order to investigate the flow field near the 
body close to and behind the base, WhithanVs 
theory requires that some form be assumed for 
the shape of the wake. Pictures of projectiles in 
flight indicate that the wake turns roughly 
through 12° at the base and then converges to 
form a cylinder of nearly constant cross section 
equal to one-half the base area. Since the exact 
form of the wake is not specified, it is of interest 
to determine the effect of various wake shapes on 
the pressure rise across the shocks in the far field. 

In order to investigate this question a shape 
obtained from equation (27) was used in which 
/l=55 and /?$//= 1/84.4, Four different wake 
shapes were then considered as shown in figure 12. 
The first is a cylinder of constant cross section 
which joins smoothly to the maximum cross 
section of the body. The two intermediate wake 
shapes separate from the body after the flow has 
turned through 15° and 30°, respectively. At 
large distances behind the body these two shapes 
become asymptotic to cylinders of constant cross 
section. The fourth wake is defined by 5^=0, 
that is, the flow is tangent to the surface every- 
where along the rear surface of the body. This 
requires the flow to turn through 38°. 

The body-shape constant was calculated to be 
the same for all four wake shapes. This result is 
believed to hold in all cases where the wake 
separates from the body behind the maximum 
section and becomes thinner thereafter but no 
general demonstration can be given. The strength 
of the tail shock, however, is affected by the shape 
of the wake as shown in reference 1. 



i 

Figure 12. — Curve showing the radius of a body shape 
with various assumed wakes. 


SEVERAL AIRPLANE AREA DISTRIBUTIONS 

References 4 and 9 indicate that the sonic-boom 
pressures associated with a nonaxisymmctric air- 
plane configuration maj^ be predicted from Whit- 
ham's theory by replacing the actual area distri- 
bution by an equivalent body of revolution 
obtained from the supersonic area rule. Whit- 
ham’s method lias recently been extended to 
lifting configurations in another manner which is 
also based on an “equivalent body” concept (ref. 
10). The five shapes shown in figure 7 arc 
equivalent-body area distributions for actual air- 
plane configurations. These shapes have a rather 
irregular area development and are included for 
comparison with the smooth shapes treated 
previously. 

Shape 2 is a Mach number 3 supersonic area 
distribution which is similar to the shapes in figure 
6. The other four shapes (1, 3, 4, and 5) are 
Macli number 1 area distributions. The bow 
shocks on shapes 1, 2, and 3 are a combination of 
three shocks which appear ahead of (x/l) maT , 
whereas the bow shocks for shapes 4 and 5 are a 
combination of two shocks. In spite of the 
diversity of the local details of these five shapes 
C b again changes very little and the pressure jump 
is largely determined by the fineness ratio. 

CONCLUSIONS 

The physical implications of the basic hypoth- 
esis for Whitman’s modification of linear theory 
(Communications on Pure and Applied Mathe- 
matics, August 1952) have been considered and 
a brief derivation given of the far-field equations 
for the bow shock. An integral w hich determines 
the effect of body shape upon the bow-shock 
overpressures in the far field has been reduced to 
a form which involves only the cross-sectional- 
area distribution and can be readily evaluated for 
body shapes for which no analytical expression is 
available. The integral has been evaluated for 
a number of families of body shapes chosen to 
investigate the effects of nose angle, fineness ratio, 
and location of maximum cross section on the 
bow r -shock overpressures. The results of these 
calculations indicate the following conclusions: 

1. In regard to body geometry, the pressure 
discontinuity in the far field is, to a first order, 
independent of body shape and depends only on 
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the fineness ratio. Local details have second- 
order effects which, in general, can be accounted 
for only direct computation of the body-shape 
constant. 

2. The calculated values of the shape constant 
which determines the effect of body geometry 
upon the pressure jump varied from 0.54 to 0.81. 
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